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Abstract 

The following paper is the continuation of our earlier considerations on cellular automata 
with Toom local rule (TCA) as the alternative to kinetic Ising systems. The arguments for 
TCA stationary states not being the equilibrium states are found in simulations. 

1 Motivation 

The Monte Carlo method applied to statistical physics denotes simulation of some stochastic 
system for which time averages restore equilibrium ensemble expectations jl], |j. Although 
the dynamical system obtained in this way is often completely artificial to the original equi- 
librium system, but much freedom in designing dynamics offers possibility to verify different 
hypothesis about micro scale interactions in the system. In this way, the Monte Carlo simula- 
tions allows to examine links between the micro dynamics and resulting equilibrium system. 
In case of the widely known Lenz model of interacting spins, the Monte Carlo method pro- 
vides the so-called kinetic Ising models. 

Why not search among other dynamical systems, such systems which also can mimic 
properties of some equilibrium system ? 

The proposition originated from cellular automata is qualitatively distinct from the men- 
tioned kinetic models. The cellular automata are complex dynamical systems. It means that 
we are given the set of local rules instead of the notion of energy and changes in the system, 
means evolution, are synchronised. The steps of synchronised update are interpreted as time 
steps. Therefore, the main goal in study such systems is to give meaning to the standard 
thermodynamic notions like energy, pressure, specific heat, temperature, etc. @, ||, |6|. In the 
case when the local rule is not reversible this goal is not obvious || [|. Especially, little is 
known about the nature of the stationary measures in the regime where there is more than 
one stationary measure. 

In the following we continue our study of cellular automata with Toom local rule (TCA) 
as the alternative to kinetic Ising systems & ||] . This model is known to exhibit non-ergodic 
properties for certain model parameters |I~0| , ||, & [jj. Therefore it can mimic the system 
undergoing the continuous phase transition. 

2 Toom cellular automata 

For every spin <7j G {— 1, +1} attached to the i node of the square lattice Z 2 we choose its 
three nearest-neighbours, named Ni,Ei,Ci as follows: 



They vote totally for the i spin state in the next time step, namely: 
Let Si = Ni + Ei + d then 

, sign E{ with probability ^(1+e) 

Ti(t + 1) = { ' (1) 

—sign Ei with probability |(1 — e) 

The parameter e G [0, 1] mimics the stochastic temperature effects: e = 1 means completely 
deterministic evolution, e = corresponds to the random rule . 

One may wish to compare Toom dynamics to the Domany rule — cellular automata with 
this rule provide the equilibrium system j|]: 

1-ef for|Ei| = 1 

sign Ej with probability < 
, ^ f 1 -e? forlEil = 3 

^ + V = \ * forlElil ^ 



-sign Ej with probability 



ef for|Ei| = 3 



and e the temperature-like parameter provides 

ef = -(1 — tanh e) and = ^(1 — tanh 3e) 

One can notice two "temperatures" playing in the Domany model. The lower temperature 
£3* is assigned to the homogeneous with respect to the spin state areas what protects them, 
while the higher temperature acts "kicking" more frequently the mixed neighbourhoods. 



3 Critical properties of TCA 

The renormalization procedure applied to equilibrium statistical mechanics systems provides 
the scaling laws, i.e. gives the relations satisfied by the critical exponents a, (3, 7, 5, u, r) ||. 
Each critical exponent characterises the power law dependence of some observable when the 
system is undergoing a continuous phase transition. Applying the standard methods for 
estimate singularities [|ll], 12] we have found the following four of the critical exponents: 

(3 — magnetisation exponent: m(e) ~ (e — e cr )^ at e > e cr 

7 — magnetic susceptibility exponent: x( £ ) ~ |e — e C r|~ 7 at e — > e cr 

v — correlation length exponent : ~ |e — e cr \~ u at e — > e cr 

rj — correlation decay at the critical point: C(i,j) ~ \i — j\~ v at e = e cr 

Our results and the methods used to obtain then are presented in Figs. 1-4. 
It is easy to check that with our estimations, namely: 

• e cr = 0.822 ± 0.002 

• = 0.12 ±0.03 

• 7 = 1.75 ± 0.03 
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• v = 0.88 ±0.02 



• 77 = 0.56 ±0.02 
the scaling lows: 

• 7 = 1/(2-7/) 

• /3 = \vr] 
are not satisfied. 

One can notice that our values obtained for and 7 are in the good agreement with those 
characterising the two-dimensional Ising system while the rest critical parameters v and r] 
are much different from the corresponding ones in the Ising system. 

Hence, TCA system seems to neither belong to the Ising class of universality 
nor be a equilibrium system . 

4 Features of TCA measures 

Statistical mechanics offers tools to investigate if a given stationary state can be represented 
by some Gibbs measure m[7|, . For a probability measure jx to be a Gibbs measure denotes 
that for any finite configuration {cja}, In /j(cta) exists and means the energy carried by the 
configuration {<ta}- The two basic features of Gibbs measures are the, so-called, quasilocality 
of interactions and the proper properties of large deviations. 
A: Quasilocality 

The idea of the quasilocality is shown in Fig. 5. Using the notation used in Fig. 5 one 
can say that nothing can arrive to some finite area A from infinity without changing the or 
configuration. In particular, if a measure \x is quasilocal then for the average magnetisation 
one have 

I < m({a A }\ ar ) > M - < m({a' A }\ ar ) > M | — > r ^oo 0. (3) 
In the following experiments we test the presence of the above property in TCA. 

I Let or be the fixed configuration evolving according to the Toom dynamics at some e. 
Let ore — the boundary configuration, be : 

a) the homogeneous configuration of all spins ±1 
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Then the (+) boundary built from the rectangle of pluses rises the probability to find 
a spin in ±1 state at the origin A = O and e' depends on the distance from the right 
and top sides of the rectangle (see Fig. 6 for details). 

Remark: If V is not of the rectangle shape then the Toom interactions enlarge V to have 
the right angle between the bottom and left sides. 
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b) the flat-interface configuration: +1 on the right and —1 on the left of the flat-interface: 



---+ + + - + + + + + 

- . ml . . + TrA - . . . . + 

- . . O . + - . m^' O . + 

- . m% . . + * - . . ..+ 
---+ + + - + + + + + 

Then the position of the flat-interface moves to the left side of the rectangle and stays 
there thanks to the periodic boundary conditions. The dependence e' on lattice site is 
the same as in the previous experiment. 

II Let ot c evolves along TCA rule at some e ou t > 0. Let the initial state for both in- 
and out- configurations is (+). Let e ou t be chosen such that it generates the stationary 
state dominated by (+) phase, e ou t ~> e cr . Then we observe: 



resulting state 

mr(l) = m(a' typical for ^+"), 
mrc(Z) = m(a' typical for fi+ ut ) 
mr(l) = m(a' typical for u £ _£ r+£ ), 
mrc(Z) = m(a' typical for n+ ut ) 
mr(l) = m(a J typical for //j_ cr ) = 0, 
mrc(i)^0 RANDOM 
mr(l) = mrc(Z) = 

These results together with the case when the initial out- phase is (— ) are presented in 
Fig. 7. Let us comment our observation as follow: 

— After adjusting the phase of the in configuration, the regions of distinct tempera- 
tures evolves independently. The two-point correlation function dies in one lattice 
unit. 

— The out-phase picks up the (+) phase from the in-system when the in-system is 
in the phase transition regime. 

— The large homogeneous structures created in the in-system undergoing the phase 
transition, propagates freely to the out configuration. Since the phase of these 
structures is random, then the phase of the out configuration changes randomly. 

— The random in-system propagates random errors outside destroying the out- phase 
magnetisation. 

Remark: The case when e out > e cr provides always mr(l) = mr<=(l) = 
Concluding: 

• If Sin < 0.76 then the conditional distribution for the magnetisation in T is independent 
of V- configuration. 

Instead, there is observed dependence in o~y c on random clusters created in o~\ if Si n is 
about the limit value and e ou t 3> e cr . 

• If E; m > 0.76 then the conditional distribution for the magnetisation depends on the 
system outside, in the sense that the probability for the spin at the origin o to take +1 



in — state out — state 
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state satisfies the inequality: 

| Prob{cr = +1 | 0T"(£m) and (Tro(e' out )} 
-Prob{ao = +1 [ cr r (e in ) and or^e^)} I > 

So that the conditional probabilities are discontinuous, what in the lattice system topol- 
ogy denotes that the property ([|) is not satisfied. The stationary measure might be 
strongly non-Gibbsian [14]. However, if Ei n > 0.84 then after thermalization time, what 



means time allowing system for adjusting the phase of the m-configuration, both areas 
evolve independently of each other. 

B. Large deviations 

For stationary measures arisen from any Markov process to be or not to be the Gibbsian 
ones is determined by the zero or non-zero value of the relative entropy density between 
different stationary measures fi and v of the system considered. li i{[i\v) > then both 
measures are non-Gibbsian |l3|, [?J. Thanks to the large deviation theorems [jl3| we have the 
powerful way to estimate i{n\v) ■ 

In case of ^_ and two stationary TCA measures corresponding to (— ) and (+) 
phases, respectively, can be extracted from the probability of the large fluctuation 

event, namely, from the probability that the large area of spins with negative magnetisation 
occurs in the stationary state described by the /i+ measure. 

From computer experiments we collect data on the magnetisation of square blocks of the 
size 1 x (. Then on the base of the formula 

ii(H-.\fi+) = lim •r='lnProb w _{m(<7j x j) < 0} (4) 

l— >oo t 



we try to estimate the limit 



i(p-\t*+) = , lim *j(m-I^+)- 

t^oo 



Fig. 8 is to present the density of the relative entropy for different lattice sizes: 

L = 60, 100, 200. Although averaging time was very large (10 000 time steps), the strong 
time dependence is noticeable. This dependence moves to the calculations of 
causing difficulties in estimating the limit I — > oo. Anyway, we tried to find out the block 
size where the relative entropy density arrives at the zero. According to the presented results 
for the lattice size L = 200 this limit should occur if the blocks are of the size I > 50. In the 
next figure, Fig. 9, we present the failure of this suggestion. What we observe in simulations 
is that with the increase of the block size the relative entropy density decay slows down. 

Moreover, we test the density of relative entropy between the stationary measure arising 
in TCA evolving little apart from the critical point, namely at £ = 0.800 and the (— ) minus 
phase, (see Fig. 10 for these results). It appears that at / = 105 the relative entropy density 
would reach the zero. But this block size is greater than the lattice size. 

Concluding, we can say that the relative entropy density between stationary mea- 
sures of Toom cellular automata evolving in the critical regime on the periodic 
lattice is positive. 

5 Conclusions 

Although the critical regime in TCA manifests itself in the way characteristic to any ther- 
modynamic system, i.e. by the rapid increase in the two-point-correlation function of spin 
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states, but the phenomena driving the system seems to be different from any equilibrium 
system. 

The TCA system undergoing the phase transition stays in the extremely dynamically 
fragile state. There is kept the sensitive balance between two processes: the processes of self- 
organising spins state to enlarge pure phase clusters and the stochastic process which destroys 
this clusters. Since the homogeneous clusters are not specially protected as it happens in the 
Domany CA, the area of cluster of one phase varies. In consequence, the correlations between 
spins are damped much stronger than in Domany CA and Ising kinetic models. However, it 
is astonishing that the relations between spin sites do not influence the critical behaviour of 
the order parameter. Therefore, one can say that the phenomenological conjecture that all 
two-dimensional ferromagnetic systems belongs to the same universality class — the class 
of Ising system [|l5| , is " partially" satisfied ( if the participation into the universality class 
could be partial). 

Testing locality of interactions we have found that at about the critical point there exists 
a boundary for the interaction to be finite. If the stochastic perturbation is sufficiently strong 
than the information from distant spins is lost. Otherwise, no matter how far away spins 
are, their influence on each other is evident. 

Our examinations on large deviation properties as well as quasilocality have hit the bound- 
ary of the finite lattice size. Therefore one can say that the system we studied does not restore 
properties of any infinite system. Such a conclusion often accompanies complex system con- 
siderations. Complex systems, existing on the, so-called, edge of order and chaos are 
known for that nobody is able to predict what kind of phenomena will arise in a system 
if one moves a little the system parameters. Studying critical properties of Toom cellular 
automata we obtain arguments for the suspicion that self-organised criticality coming from 
the dynamic equilibrium provides systems qualitatively distinct from the thermodynamic 
systems [[l?]]. 
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Figure 1: Estimates for (3 in Toom CA. Data collected in 430, 280, 120 experiments with lattices: 
L = 100, 150, 200 respectively. 

Figure 2: Estimates for 7 in TCA from the relation k =< m 2 > — < m > 2 . Data collected in 
430, 180, 120 experiments with lattices: L = 100, 120, 150 respectively. 

Figure 3: Estimates for \jv in TCA on the base of finite size theory from the maximal slope 
of derivative of the 4th order magnetisation cumulant as well as from the maximum slope of the 
logarithm derivatives of: absolute value of magnetisation, square magnetisation, fourth power 
of magnetisation. Data collected in 1400, 860, 660, 430, 280, 120 experiments with lattices: 
L = 20, 30, 60, 100, 150, 200 respectively. 

Figure 4: The decay rate n for the two-point correlation function of magnetisation obtained on 
the lattice with L = 100 at different e . 

Figure 5: The illustrative definition of quasilocality of interaction in the lattice system. A and V 
are any finite subsets of a lattice C . The idea is that in case of large V there is not observable 
influence on the configuration in A coming from from any configuration a T c which is outside to 
the fixed configuration of a r . 

Figure 6: TCA with (+) boundary added after reaching stabilisation. Distinct curves correspond 
to different values of stochastic perturbation e. With e' we measure the extra magnetisation 
observed at the origin O in stationary TCA. 

Figure 7: Interaction between two TCA systems: in and out distinct from each other because of 
different stochastic parameters: e in and e out , suitable, a) in initial state is (+)- phase, e in = 0.86, 
out initial states are changed. Labels in the figure correspond to different values of e out . b) 
Propagation of self-created clusters of one phase from in area into the out configuration observed 
as random changes in magnetisation of out state. 

Figure 8: The relative entropy density between a + and yU_ stationary measures in the critical 
regime of TCA versus block size for lattices L = 60, 100, 20. 

Figure 9: The relative entropy density between u + and /z_ stationary measures for large blocks. 

Figure 10: Density of the relative entropy between stationary measures /i_ and v — the measure 
for the system moved a little from the critical point. 
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